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Contrasting  O/X-mode  heater  effects  on  O-mode  sounding  echo  and 
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[i]  The  effects  on  the  ionosphere  of  powerful  O-mode  and 
X-mode  HF  pump  waves,  modulated  3  minutes  on  and 
1  minute  off,  were  explored.  The  experiments  were 
monitored  using  the  digisonde  and  magnetometer 
located  at  the  HA  ARP  facility.  The  results  show  that  the 
virtual  heights  of  the  O-mode  sounding  echoes  shifted 
down/up  as  the  0/X  mode  heater  was  turned  on;  the 
ionosphere  also  moved  downward/upward  accordingly. 
Enhanced  spread-F  was  also  observed  in  O-mode  heater-on 
periods.  Heater-induced  magnetic  pulsation  was  observed. 
Its  intensity  increased  progressively  in  the  heater  on/off 
sequence  and  X-mode  heater  was  more  effective  than 
O-mode  heater  in  the  generation  of  magnetic  pulsation.  In 
the  last  X-mode  heater-on  period,  when  the  magnetic 
pulsation  reached  the  highest  level,  pc  3  pulsations,  with 
increasing  intensity  were  also  observed.  Citation:  Kuo,  S., 
W.-T.  Cheng,  A.  Snyder,  P.  Kossey,  and  J.  Battis  (2010), 
Contrasting  O/X-mode  heater  effects  on  O-mode  sounding  echo 
and  the  generation  of  magnetic  pulsations,  Geophys.  Res.  Lett.y 
37,  LOllOl,  doi:10.1029/2009GL041471. 

1.  Introduction 

[2]  A  major  facility  for  conducting  ionospheric  heating 
experiments  is  available  in  Gakona,  Alaska,  as  part  of  the 
High  Frequency  Active  Auroral  Research  Program 
(HA  ARP)  [Kossey  et  al.^  1999].  HA  ARP  HF  transmitting 
system  consists  of  1 80  crossed  dipole  antennas  arranged  as 
a  rectangular,  planar  array.  Each  antenna  element  radiates 
circularly  polarized  wave  up  to  20  kW  in  the  frequency 
band  from  2.8  MHz  to  1 0  MHz,  over  which  the  antenna  gain 
increases  from  15  dB  to  30  dB  with  increasing  radiating 
frequency.  It  is  anticipated  that  HF  heating  of  the  ionosphere 
will  make  a  significant  impact  by  both  collision  and  anom¬ 
alous  heating  processes. 

[3]  In  heating  experiments,  the  HF  heater  wave  is  nor¬ 

mally  transmitted  in  the  X-mode  (left-hand  (LH)  circular 
polarization),  or  in  the  O-mode  (right-hand  (RH)  circular 
polarization),  with  the  heater  frequency  ujq  <  the 

plasma  frequency  of  the  F  peak.  These  heating  waves 
impact  the  ionosphere  through  different  processes.  One 
process  is  to  excite  parametric  instabilities.  However,  this 
requires  that  the  heating  wave  be  accessible  to  the  spatial 
regions  where  the  parametric  coupling  conditions  can  be 
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matched  [Fejer^  1979;  Stenflo,  1985].  Consequently,  the 
region  directly  impacted  by  instabilities  is  limited  to  close 
to  the  reflection  layer  of  the  heater.  In  the  present  work, 
an  experiment  was  conducted  to  explore  the  effects  of 
O/X-mode  heating  on  a  broad  region  of  the  ionosphere  by 
examining  digisonde  and  magnetometer  records.  ULF 
waves  can  be  produced  by  the  HF  heater  via  fi  lamentation 
instability  [Kuo  and  Lee,  1983]  or  other  dynamo  processes 
[Papadopoulos  and  Change  1985].  Filamentation  instabil¬ 
ity  [Kuo  and  Schmidty  1983]  simultaneously  generates 
large-scale  field-aligned  density  irregularities  (FAIs)  as 
has  been  observed  in  a  recent  HAARP  heating  experiment 
[Kuo  et  al.y  2009].  Hence,  the  magnetometer  records  may 
be  used  to  infer  the  presence  of  heater-induced  large-scale 
FAIs  in  addition  to  the  scintillation  method  used  by  others 
[Frey  et  al.^  1984;  Secan  et  al.,  1997]. 

2,  Experimental  Observations 

[4]  On  March  18,  2009  from  UTC  20:45  to  UTC  21:08, 
from  UTC  21:35  to  UTC  22:08,  and  from  UTC  22:11  to 
UTC  22:50,  using  the  HAARP  transmitter  facility  at 
Gakona,  AK,  at  full  power  (3.6  MW)  experiments  were 
conducted  with  HF  heater  waves  of  3.2  MHz  directed  along 
the  geomagnetic  zenith.  From  20:45  to  21:08,  O-mode  and 
X-mode  heating  modulated  with  3  minutes  on  and  1  minute 
off  were  transmitted  alternatively.  The  rectangular  modula¬ 
tion  wave  had  6  cycles  in  the  24  minutes  experimental 
period.  In  the  next  two  time  periods,  only  X-mode  pump 
waves  were  transmitted.  From  21:35  to  22:08,  the  heater 
was  modulated  with  3  minutes  on  and  3  minutes  off  for 
6  cycles,  and  then  the  modulation  was  changed  to  3  minutes 
on  and  1  minute  off  for  10  cycles  in  the  period  from  22:1 1 
to  22:50. 

[5]  The  experiments  were  monitored  by  the  HAARP 
digisonde  and  magnetometer.  Because  the  HF  heating  wave 
is  in  the  frequency  band  of  the  digisonde,  it  can  be  recorded 
by  the  digisonde  receiver  as  a  false  return  with  much 
stronger  intensity  that  degrades  the  image  resolution  of  the 
true  returns.  To  avoid  this  occurence,  the  ionogram  is 
blanked  in  a  frequency  band  around  the  heater  frequency. 
Specifically  for  these  examples,  the  returns  in  the  band 
approximately  from  3  to  3.4  will  be  blank.  The  exact 
ionogram  bandwidth  around  the  heater  frequency  is  deter¬ 
mined  by  a  combination  of  inonsonde  functions  such  as 
automatic  gain  control,  signal  threshold  determination, 
and  the  signal  detection  threshold  of  the  receiver. 

2.1.  Heater  Induced  Virtual  Height  Variation 

[6]  The  ionograms  in  the  frequency  bands  of  2.5  to 
2.9  MHz  and  3.55  to  3.95  MHz  were  then  examined  to 
extract  evidence  of  generating  large-scale  FAIs  by  the  HF 
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Figure  1.  Sequence  of  ionograms  of  0-mode  sounding  echoes  in  the  frequency  band  of  2.52  to  2.73  MHz,  recorded  from 
20:47  to  21:08  with  3.2  MHz  0/X-mode  HF  pump  waves,  modulated  with  3  minutes  on  and  1  minute  off,  transmitted 
alternatively. 


heater.  The  digisonde  was  run  in  a  relatively  fast  mode 
that  produced  ionograms  during  the  heater  on  and  off 
periods  for  purposes  of  comparison. 

[7]  We  first  examined  the  ionograms  recorded  in  the  first 
experimental  period  from  20:45  to  21:08.  The  effect  of  the 
HF  heating  on  the  virtual  heights  of  the  0-mode  sounding 
echoes  in  the  frequency  band  of  2.52  to  2.73  MHz  was 
explored  by  directly  comparing  the  ionograms.  The  results 
of  six  consecutive  heater-on/off  sequences,  presented  in 
Figure  1,  show  that  the  virtual  heights  of  the  echo  shift  down 
(or  up)  as  the  O  (or  X)  mode  heating  is  turned  on.  The 
recording  of  the  ionogram  was  started  at  the  beginning  of 


a  minute,  and  the  time  span  for  a  complete  scan  from  1  to 
6  MHz  was  about  37  seconds.  Thus,  echoes  in  the  2.52  to 
2.73  MHz  band  of  the  off  period  ionograms  were  recorded  at 
about  1 1  to  1 3  seconds  after  the  heater  was  turned  off.  As  seen 
in  Figure  1,  although  the  virtual  heights  of  the  echoes  in  each 
off  period  ionogram,  are  not  much  different  from  those  in 
the  preceding  on  period  ionogram,  which  suggests  that  the 
ionospheric  recovery  time  is  longer  than  13  seconds,  the 
change  of  the  color  code  of  the  echoes  has  a  significant 
meaning.  In  the  plots,  red  and  blue  echoes  represent  negative 
and  positive  Doppler  shifted  vertical  0-mode  echoes,  i.e., 
echoes  from  up-moving  and  down-moving  ionosphere.  In 


21:35  21:38  21:47 


21.-50  2159  22.D2 


Figure  2.  Sequence  of  ionograms  of  0-mode  sounding  echoes  in  the  frequency  band  of  2.68  to  2.87  MHz,  recorded  from 
21:35  to  22:02  with  3.2  MHz  X-mode  HF  pump  wave,  modulated  with  3  minutes  on  and  3  minutes  off. 
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Figure  3.  Effect  of  0-mode  HF  pump  wave  on  the 
(a)  virtual  height  and  (b)  intensity  of  0-mode  sounding 
echo. 


one  cycle  with  0-mode  heater  on/off  and  then  X-mode  heater 
on/off,  the  color  change  of  the  majority  of  the  echoes,  from 
blue  to  red  and  then  red  to  blue,  indicates  that  0-mode  heater 
moves  the  ionosphere  down  and  X-mode  heater  moves  it  up; 
and  in  the  off  period,  ionosphere  moves  up  (or  down)  to 
recover  the  prior  down  (or  up)  shift  by  the  O  (or  X)  mode 
heating.  Figure  1  also  shows  that  0-mode  heater  caused  an 
artificial  spread- F,  which  persisted  even  after  the  heater  was 
turned  off  for  more  than  the  1 3  seconds  corresponding  to  the 
time  delay  in  transmitting  sounding  pulses  in  the  frequency 
band  of  2.52  to  2.73  MHz. 

[s]  The  fact  that  X-mode  heater  moves  the  ionosphere  up 
can  be  seen  clearly  during  the  second  experimental  period 
from  21:35  to  22:08.  Three  sets  of  on/off  ionograms  in  the 
frequency  band  of  2.68  to  2.87  MHz  are  presented  in 
Figure  2  to  show  the  consistent  effect  of  the  X-mode  heater 
on  lifting  the  ionosphere. 

[9]  This  contrasting  effect  of  O  and  X  mode  heating  on 
the  virtual  heights  of  the  0-mode  sounding  echoes  also 
appears  in  the  frequency  band  of  3.55  to  3.95  MHz  of  the 
ionograms.  This  is  exemplified  by  showing  two  cases 
analyzed  from  two  sets  of  heater  on/off  ionograms,  which 
were  recorded  in  20:03/20:04  and  21:35/21:38.  The  fre¬ 
quency  dependencies  of  the  virtual  heights  of  the  highest 
amplitude  0-mode  sounding  echoes  in  the  two  cases, 
corresponding  to  O  and  X-mode  heater  on/off,  are  plotted 
in  Figures  3a  and  4a,  respectively.  As  shown,  the  virtual 
heights  of  the  echoes  were  also  downshifted  by  the  0-mode 
heater  (Figure  3a)  and  upshifted  by  the  X-mode  heater 
(Figure  4a).  This  contrasting  trend  is  also  identified  visually 
from  the  corresponding  ionograms  presented  in  Figures  3b 
and  4b,  in  which  the  point  of  the  intensity  peak  at  each 
sounding  frequency  is  marked  by  a  black  dot.  It  is  noted  that 
every  echo  in  the  ionograms  of  Figures  1  -4  is  plotted  with 
the  same  size  rectangle  regardless  of  its  amplitude. 

2.2.  Excitation  of  Geomagnetic  Micropulsation 

[10]  The  induced  magnetic  field  variation  was  monitored 
by  the  fluxgate  magnetometer  located  at  Gakona,  AK.  The  1 
sec  resolution  data  from  UTC  20:45  to  UTC  21:08  were 
averaged  for  over  1 0  sec  periods  providing  1 8  average  values 
in  each  heater  on  period  (180  sec)  and  6  average  values  in 
each  heater  off  period  (60  sec).  These  average  values  contain 
both  the  background  magnetic  field  as  well  as  any  variation 
caused  by  changes  of  the  electrojet.  A  four  point  algorithm  is 


developed  to  extract  the  heater-induced  component  from  the 
average  data.  This  algorithm  extrapolates  heater- independent 
magnetic  field  in  a  heater-on  period  by  using  the  heater-off 
values  at  the  two  ends  of  the  heater-on  period.  Each  point  is 
extrapolated  by  the  three  adjacent  points,  one  in  the  front  and 
two  in  the  back.  The  front  point  is  unknown,  thus  the 
extrapolation  is  done  from  both  directions  simultaneously. 
Consider  one  on  period  and  its  two  adjacent  off  periods; 
taking  only  2  average  values  in  each  heater  off  period  and 
label  the  22  average  values  (2  +  1 8  +  2)  in  these  three  periods 

(off-on-ofT)  by  Aj,  j  =  1,2 . ,  22,  where  (Ai,  A2)  and 

(A21,  A22)  are  the  average  values  in  the  front  and  back  off 

periods.  Let  Aj  =  Af  +  ABj,  where  ABj,  j  =1,2 . 22,  are 

the  average  magnetic  field  produced  by  the  heater  in  each 
10  sec  slot,  where  (ABi,  AB2)  and  (AB21,  AB22)  are  zero. 
Using  the  relations,  deduced  approximately  from  the 
Runge-Kutta  algorithm  [Conte  and  de  Boor^  1972]  by 
representing  the  slope  functions  in  the  algorithm  in  terms 
of  finite  differences  between  points,  An  —  An+I  —  y3An+2  + 
ysAn-i  ^  An-1  -  y3An-2  +  l^An+l,  yield 

AB„  =  A„-(I/80){[n(22-n)-f  (n-  I)(2I  -  n)]A2 

-[(n-2)(22-n)  +  (n-3)(2l-n)]Ai} 

-  (I/80){[(n  -  I)(23  -  n)  +  (n  -  2)(22  -  n)lA2, 

-  [(n  -  I)(2I  -  n)  +  (n  -  2)(20  -  n)]A22} 

forn  =  3,4 . 19,20;  (I) 


[11]  We  then  use  (1)  to  process  the  averaged  magnetom¬ 
eter  data.  The  results  reveal  the  evolution  of  the  10  sec 
average  magnetic  field  induced  by  the  heater  in  one  on- 
period. 

[12]  The  results  of  ABh  (pink  curve),  ABq  (blue),  and 
ABz  (green)  in  three  directions  (H,  D,  and  Z),  from  UTC 
20:45  to  21:09,  are  presented  in  Figure  5.  As  shown,  the 
magnetic  field  fluctuation  induced  by  the  X-mode  heater  is 
larger  than  that  induced  by  the  0-mode  heater;  the  intensity 
and  oscillatory  features  of  the  induced  magnetic  fields,  in 
both  heating-mode  cases,  develop  with  the  duration  of  the 
heating  cycles.  The  peak  intensity  of  ABp  reaches  1.5  nT 
and  2  nT  in  the  last  0-mode  (21:01  to  04)  and  X-mode 
(21:05  to  08)  heating  cycle.  In  the  ABh  and  ABz  compo¬ 
nents,  the  development  of  oscillatory  feature  is  more  pro¬ 
nounced.  Furthermore,  ABh  in  the  last  X-mode  heating 
cycle  has  evolved  into  pc  3  pulsations  with  increasing  inten- 


Figure  4.  Effect  of  X-mode  HF  pump  wave  on  the  (a)  vir¬ 
tual  height  and  (b)  intensity  of  0-mode  sounding  echo. 
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Figure  5.  H  (pink),  D  (blue),  and  Z  (green)  components 
of  the  magnetic  field  fluctuation  induced  by  the  HF  pump 
wave. 

sity.  The  polarities  of  ABh  and  AB^  in  X-mode  heating 
cycles  are  basically  opposite  to  each  other;  on  the  other  hand, 
ABh  ABd  have  the  same  polarity  in  0-mode  heating 
cycles.  These  contrasting  differences  in  the  induced  magnetic 
fluctuations  suggest  that  the  density  irregularities  generated 
by  the  O  and  X-mode  heaters  are  quite  different  in  the  char¬ 
acteristic  features  including  polarization,  scale  size,  and 
magnitude. 

3.  Discussion 

[13]  The  0-mode  heater  is  reflected  at  a  height  near  the 
electron  plasma  resonance  layer,  where  the  electron  plasma 
frequency  o^po  =  and  is  located  higher  than  the  upper 
hybrid  resonance  layer  a^pu  =  —  0^)*^^.  Therefore, 

0-mode  heating  wave  is  accessible  to  the  spatial  regions 
where  many  parametric  coupling  conditions  can  be  matched 
[Fejefy  1979;  StenJIo,  1985].  Parametric  instabilities  convert 
HF  heating  wave  to  plasma  waves  which  heat  the  ionosphere 
locally  near  the  HF  reflection  height.  Such  anomalous  heating 
forms  a  local  heat  source  which  transmits  heat  downward 
through  heat  conduction  along  the  magnetic  field.  Conse¬ 
quently,  the  ionosphere  moves  downward  and  a  thermal 
instability  [Kuo  and  Djuth,  1 988]  is  also  excited  to  generate 
large-scale  FAIs  which  enhance  spread-F  as  seen  in  Figure  1 . 

[  14]  On  the  other  hand,  the  X-mode  heater  will  be  reflected 
at  a  height  having  a  plasma  fi^equency  of  -  [cjo(^o  ” 
Oy)]*^  that  is  below  the  electron  plasma  resonance  layer  as 
well  as  below  the  upper  hybrid  resonance  layer.  Although  the 
X-mode  heater  cannot  effectively  introduce  a  similar  anom¬ 
alous  heating  in  the  F-region  it  heats  the  ionosphere  through 
the  collision  process  more  effectively  than  does  the  0-mode 
heating.  The  quiver  speeds  |Vqi|  of  electrons  in  the  0/X- 
mode  heating  wave  electric  field  E  are  given  by  |Vq±|  = 
|eE/m|[(a;o  ±  Og)^  +  ,  where  Og  and  u  are  the  electron 

cyclotron  and  collision  fi^equencies.  Thus  the  ratio  ^  of  the  X 
to  0-mode  collision  heating  rates  is  obtained  to  be  ^  =  [(u;o  + 
fief  +  +  A  ^  (Wo  +  ne)^/(Wo  -  »  1  • 

The  collision  heating  process  is  non-local  and  moves  the 
ionosphere  upward. 

[15]  The  electron  thermal  pressure  force,  arising  from  the 
differential  ohmic  heating  on  electrons  [Gurevich,  1976; 
Fejer,  1 979]  by  the  combined  wave  electric  field  of  the  HF 
heater  and  the  excited  high  frequency  sidebands,  pushes 
electrons  to  form  laige-scale  FAIs,  which  in  turn  break  up 


the  HF  pump  wave  by  filamentation  instability  [Kuo  and 
Schmidt,  1983].  The  electron  thermal  pressure  force  Fn  also 
gives  rise  to  an  Fh  x  Bq  drift  motion  in  the  electron  fluid 
and  induces  a  net  electron  drift  current  flowing  perpendicular 
to  both  the  background  magnetic  field  Bq  and  the  spatial 
variation  direction  of  the  produced  density  irregularities. 
Therefore,  magnetic  field  pulsations  are  excited  simulta¬ 
neously  with  large  scale  density  irregularities  by  the  filamen¬ 
tation  instability  [Kuo  and  Lee,  1 983;  Lee  and  Kuo,  1 985]. 
Pc  5  pulsations  have  been  observed  in  Tromso  heating 
experiments  [Stubbe  and  Kopka,  1981]  and  ULF  waves 
have  been  observed  by  Chang  et  al.  [2008]  in  HAARP 
heating  experiments.  However,  it  is  noted  that  the  present 
experiment  was  conducted  around  local  solar  noon  during  a 
geomagnetic  quiet  time,  thus  the  generation  of  magnetic 
pulsations  is  electrojet  independent.  The  collision  heating 
by  the  X-mode  heating  is  more  effectively  than  that  by  the 
0-mode  heating  and  it  explains  the  observation  in  Figure  5 
that  the  heater-induced  magnetic  pulsations  are  more  intense 
in  X-mode  heater-on  periods.  The  heater  was  modulated  with 
3  minutes  on  and  1  minute  off  so  the  generation  of  Pc  3 
pulsations  in  the  X-mode  hcater-on  period  was  not  directly 
related  to  the  modulation  period.  Pc  3  pulsations  could  be 
generated  by  the  current  pulsations  in  closing  the  loops  of 
currents  driven  by  heater-induced  thermal  pressure  force. 
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